Rotational states in deformed nuclei: An analytic approach 
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The consequences of the spontaneous breaking of rotational symmetry are investigated in a field 
theory model for deformed nuclei, based on simple separable interactions. The crucial role of the 
Ward-Takahashi identities to describe the rotational states is emphasized. We show explicitly how 
the rotor picture emerges from the isoscalar Goldstone modes, and how the two-rotor model emerges 
from the isovector scissors modes. As an application of the formalism, we discuss the Ml sum rules 
in deformed nuclei, and make connection to empirical information. 
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I. INTRODUCTION 

The common approach to describe deformed nuclei is 
based on the mean field (Hartree) approximation, where 
the rotational symmetry of the Hamiltonian is sponta- 
neously broken [H 0]. As a consequence of the broken 
symmetry, Goldstone poles emerge in the Bethe-Salpeter 
(BS) equation (or, equivalently, the RPA equation) for 
a particle- hole pair[3|, similar to the case of infinite 
systems Q. For finite systems with axial symmetry, these 
intrinsic zero modes lead to the picture of collective rota- 
tion of the whole system around an axis perpendicular to 
the symmetry axis, thereby restoring the symmetry of the 
original Hamiltonian The corresponding rotational 
band is characterized by finite excitation energies. The 
self consistency relations for the deformed mean fields [6] 
provide the necessary conditions for the existence of the 
Goldstone poles and the ground state rotational band. 
These features, which will be elucidated in this paper by 
using a simple field theory model for nucleons, were dis- 
cussed transparently in models based on bosonic degrees 
of freedom Q , and form the basis for the construction of 
effective theories for deformed nuclei Q. 

Besides these collective rotations of the whole sys- 
tem, which can be called isoscalar rotations, there ex- 
ist also rotational modes of isovector character, i.e., ro- 
tational vibrations of protons against neutrons [H, Hoj . 
These so-called scissors modes, which were predicted 
originally in the 2-rotor model II l[ and further investi- 
gated by using sum rule methods [12| and the Interact- 
ing Boson Model [l^. can be excited by the isovector or- 
bital part of the Ml operator, and decouple automati- 
cally from the isoscalar rotational modes if the self con- 
sistency relations are satisfied [fj, Recent experimen- 
tal investigations [TBI on the scissors modes have con- 
centrated on magnetic sum rules, which are very impor- 
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tant because they provide the connection of the collective 
modes to quantities like the effective orbital g-factors or 
moments of inertia, which can be determined by other 
independent analyzes (l7j. The experimental and theo- 
retical works on the scissors modes and other magnetic 
dipole modes at higher excitations energies [18j are sum- 
marized in a recent extensive review (l9|. 

The purpose of the present paper is threefold: First, we 
wish to elucidate the importance of the Ward-Takahashi 
identities [20j to describe the intrinsic isoscalar Goldstone 
modes and the ground state rotational band. In particu- 
lar, we wish to show that the results derived for example 
in Ref.@ can be obtained rather elegantly by using the 
Ward-Takahashi identities without explicit reference to 
single particle wave functions or truncations of the model 
space. Second, we wish to discuss how the isovector scis- 
sors mode, which corresponds to a RPA solution with 
finite energy, leads to the two-rotor picture in a simple 
model calculation. (To the best of our knowledge, a sim- 
ple and direct derivation of the rotor picture from the 
intrinsic Goldstone modes, and of the two-rotor picture 
from the isovector scissors modes, has not yet been pre- 
sented.) Third, we wish to investigate the inverse energy 
weighted and energy weighted Ml sum rules, and discuss 
the relation to recent experimental works on the scissors 
modes. For these purposes, we will use a simple field the- 
ory model based on a separable quadrupole-quadrupole 
(QQ) interaction in the BS (RPA) framework, and con- 
sider only the orbital motion of the nucleons. Our main 
interest here is the physics of the rotational modes at low 
energies in well deformed heavy nuclei, and there is evi- 
dence both experimentally [l9| and theoretically [l(| that 
these low energy modes are basically of orbital character. 

It should be noted here that separable interactions 
have widely been used to investigate deformed nuclei in 
the RPA|21ri23|. In order to perform quantitative calcu- 
lations, it is well known that pairing plays an important 
role. The Nambu-Gorkov formalism [24j actually pro- 
vides a simple way to incorporate the effects of pairing 
into the properties of quasiparticles. Our purpose here, 
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however, is to gain analytic insights into the physics be- 
hind the rotational states and the associated sum rules 
in the most simple and transparent way. Therefore, the 
pairing effects will not be included in the formulas, and 
numerical results will not be presented in this paper. Our 
hope is that the analytic approach described here can be 
extended to more general types of interactions. 

In Sect. 2 we will formulate the mean field approx- 
imation and the BS equation for our model. In Sect. 
3 we will use the Ward-Takahashi identities for angular 
momentum conservation to derive several important rela- 
tions and low-energy theorems for Green functions. The 
connections to observables will be established in Sect. 4, 
where we will discuss basic properties of transition matrix 
elements, and in Sect. 5, where the Ml sum rules will be 
derived. In Sect. 6 we will discuss the derivation of the 
rotor picture from the isoscalar Goldstone modes, and 
of the 2-rotor picture from the isovector scissors modes. 
For the derivation of the ground state rotational band, 
no further approximations are necessary, but in order to 
derive the two-rotor picture we still find it necessary to 
assume a harmonic oscillator potential for the spherical 
part of the mean field. A summary and an outlook are 
given in Sect. 7. 

II. THE MODEL 

The Hamiltonian of the model which we will use in this 
paper is given by 

H = h 0p + h 0n + ^ (Qt • Q p + Qt . Q n ) 
+ *f-(Qt'Qn + Qn • Qp) • 

Here h 0T = / d^t (x) H ( X ) Mx ), where 

Hq(x) = — A/2M + Uo(r ) with M the nucleon mass 
and Uq(t) some spherical mean field. The quadrupole 
operator is defined by 

Q 1 ; = J d 3 xri(x)Q K (x)Mx) (T=P,n), (11.2) 

where Q K (x) — r 2 YiK(i) 1 , and the products ■ Q 
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in (|ILT|) are denned as Q ] ■ Q = ^ Q K ' I Q K . For the 

K=-2 

coupling constants of the QQ force we assume \nn — Xpp 
and Xpn = Xnp- Because we treat protons and neutrons 
as separate particles, the interaction in 1|) is a mixture 
of pure isoscalar (xpn = Xpp) an d pure isovector (xpn = 
—Xpp) type interactions. 



A. Mean field approximation 

The mean field approximation is formulated as usual 
by adding and subtracting a term —J2 T PtQ t + C in 
(III. II) . where the parameters f3 T and the constant C will 
be determined later by the requirement of self consis- 
tency. In order to avoid mathematical ambiguities in the 
low energy theorems to be discussed later, we also add the 
terms — ^2 t £tQ t , which explicitly break the rotational 
symmetry. (In the final results the symmetry breaking 
parameters e T will be set to zero 2 .) In this way we obtain 

H = H 0p + H 0n + C + {(f3 p - e p ) Q° p + (/3„ - e„) Q° n 

+ *f (Ql ■ Q P + Ql ■ Qn) 

+ ^f- (Ql- Qn + Qt- Q P ) - c], (n.3) 

where 

Hot = J d 3 xV4(x) {-^ + U (x) - l3 T Q°(x)^j i> T (x). 

(II.4) 

We now assume that the rotational symmetry is spon- 
taneously broken, i.e., that only the K = component 
of the quadrupole operator has a finite ground state ex- 
pectation value: 

Q° T = (Q° T )+ : Q° T : , (II.5) 

where the dots in the second term denote normal order- 
ing. We require that the part {...} in (|II.3|) becomes 
a "true" residual interaction, i.e; when (|II.5[) is inserted 
into (|II.3p . this part has neither terms linear in : nor 
constant (c-number) terms. The first requirement leads 
to the self consistency relations 

ftp — £p ~ XppiQp) ~ Xpn (Q 

Pn = e n -Xnn{Q n )-Xn P {Q° p ), (H.6) 

and the second requirement determines C as 

C = -*f {{Ql? + {Qn) 2 ) - Xpn(Q Q p )(Q° n ). (H.7) 

Self consistency implies that the expectation values 
inEqs. (|II.6|) themselves depend on Bt. The Hamiltonian 
finally becomes 

H = H 0p + H 0n + C + J2^Y ( : Q ^ : ) ' Qp • 

Tp 

(11.8) 



1 For clarity, we will also use the notations Q^(x) or (x) to 
indicate whether the quadrupole field refers to protons or neu- 
trons. In this paper, the labels t, p, A stand for protons (p) or 
neutrons (n). If a sum over those labels is involved, it will be 
indicated explicitly. 



2 These symmetry breaking parameters e should not be confused 
with the single particle energies, for which we will use the symbol 
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FIG. 1: Graphical representation of the particle-hole interac- 
tion kernel Eq. (|II.9[) . In this and all following diagrams, time 
can be visualized to run from left to right. 




FIG. 2: Graphical representation of the particle-hole T- 
matrix Eq. (IH.10|l in the ladder approximation. 



If one employs a harmonic oscillator potential Uq (r) = 
(Afw 2 /2) r 2 , it is sometimes convenient to define dimen- 
sionless deformation parameters f3 T by f3 T = AIuj 2 f3 T , be- 
cause then one can express the sum Uq(x) — /3 T Q°(x) in 
(111.41) as a deformed harmonic oscillator potential and ap- 
ply standard methods of the Nilsson model[25j]. Except 
for the last parts of this paper (Sect. 6.B), we will keep 
the discussions general without specifying the form of Uq. 



B. Collective excitations 

Here we consider the Bethe-Salpeter (BS) equation, 
which is equivalent to the RPA equation, for a particle 
and a hole in the K = 1 channel. 3 

The residual interaction in pi.8[) is separable in coor- 
dinate space, see Fig.l. The corresponding Feynman rule 
for the particle-hole interaction kernel is given by 



K rp (x',x) = -iQKx^XrpQlHx). 



(II.9) 



The inhomogeneous BS equation then reads (see Fig. 2) 
T Tp (x',x;u) = K Tp (x\x) + J dV| dhj / 1^ £ 

K T \ (x\ y') S\ (y, y; k a + u) S\ (y, y'; k ) T Xp (y, x; oj) . 

(11.10) 

Here we work with a mixed representation of the Feyn- 
man propagator: 

a l , s _ \ - <j>r a {x')(l>\ a {x) \ - (jjTijx'Wl^x) 

O t {X , X, UJ) — > — — h / — , 

*■ — ' u! — e Ta + id L — ' uj — e T i — id 

= S TP + S TH , (11.11) 

where 4> Ta {x) and e Ta are the eigenfunctions and eigen- 
values of Hq t for particle (P) states, and c/) T i(x), e T i de- 
note the corresponding quantities for hole (H) states 4 . 
Inserting the kernel Eq. (|II.9|l and the ansatz 

T tp (x',x;lu) = -iQlix'^^Q^ix) (11.12) 



into pi,10[) . we obtain the following simple matrix equa- 
tion for the reduced T-matrix: 



*M = X-%T(w)t(u)), 
. , 1 1 

H w ) = i I ^ X = X 



(11.13) 



1 + XttM l + 7r(o;)x 

Here the matrices in charge space have the form 

f = (tpp tpn \ _ ( Xpp Xpn \ n = f^P 

\t n p t nn J yXnp Xnn J y ^"n 

(11.14) 

and the proton and neutron "bubble graphs" tt t (uj) are 
given by (see Fig. 3 and Appendix A) 

x [Q 1 T \x')S T (x',x;k +uj)Q 1 t (x)S t (x,x';k )] , 

(11.15) 



-2 J2 IWQ'M 



(az)67 



Here oj a i — e a — 6j are the non-interacting particle-hole 
energies. 

From (|II.13j) . the poles of the T-matrix (oj 2 = w 2 ) are 
determined by the equation 

Det(l + X7r(cj))) 

= (1 + ir p (oj)xpp) (1 + 7T n (w)Xnn) - Xp„7i>M7r™(w), 
= 0. (11.16) 



3 The K = — 1 channel is degenerate with the A" = 1 channel, 
while the K = and the = ±2 channels have different ener- 
gies. There is no mixing of those channels for the case of axial 
symmetry. We consider the collective K = ±1 states here, be- 
cause they correspond to the rotational states which are of main 
interest in this paper. 



4 Notations like a £ r (or i £ r) indicate that the single-particle 
state a (or the single-hole state i) is a proton (r = p) or neutron 
(r = n) state. We also remark that the states i in 1)11. llj l are 
actually the time-reversed of the occupied single particle states. 
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FIG. 3: Graphical representation of the bubble graph 
Eq. (HLT51) . 



FIG. 5: Graphical representation of the homogeneous BS 
equation, Eq. ([II.24j) . The double line represents a collective 
state. 



xp ' 



CO 2 - co^ 



Here the prime indicates differentiation w.r.t. u> 2 , i.e., 



<(u,) = = 2 J] IHQrl 



dw 



(az)g7 



2 ' 



(11.23) 



FIG. 4: Graphical representation of the pole behavior of the 
particle- hole T-matrix, Eq. (|H.18|) . The double line represents 
a collective state. 



It is straight forward to use Eq. (|II.13j) to determine 
the pole behavior of the T-matrix. The results for the 
reduced and the full T-matrix are (see Fig. 4 for the full 
T-matrix) 



t Tp (oj) 
T tp (x',x;uj) 



u 2 ^l N T (uj n )N p {u) n ) 
us 2 — lo 2 + iS 

up- — ui 2 + id 



Here the vertex functions for the collective K 
n (excitation energy u> n ) are given by 

T?(x)=Qy(x)N T (u n ), 

with the normalization factors determined from 5 

N n (ui n ) 2 



= 1 



N p (co n ) 2 py H ' ny n > N p (co n ) 2 



N n {u) n ) _ 1 + XppKpjUn) _ -XpnKp{LL>„) 

K) i + K) 



11.17) 
11.18) 

state 

11.19) 
11.20) 

11.21) 
11.22) 



It is also straight forward to derive the above forms of the 
vertex functions (except for the overall normalization) 
from the homogeneous BS equation: Inserting the pole 
behavior (|II.18|) into Eq. pi.lO[> and taking the limit 
uj 2 — > uj^, one obtains the homogeneous BS equation (see 
Fig. 5): 



r?(a:) 



x Q\[x) Xr P Q 1 p \y)S P {y 1 y'; k + u)S p (y', y; k ) T n p (y'). 

(11.24) 



Inserting here the ansatz T™(x) = Q\ (x)N T (w), one 
obtains the following matrix equation for the normaliza- 
tion factors: 



N n (u)J 



XppKpiuj) XpnK n (u)\ ( N p (ui) 

Xnp^p (uj) J \N n ((jj) 



(11.25) 



This equation again leads to the pole equation (|II.16[) 
with solutions uj = uj n , and to the relation (|II.22|) . 



5 We normalize the vertex functions T as the residues at the poles 
in up, which corresponds to "covariant normalization" in rela- 
tivistic field theory. The overall sign is chosen so that for a pure 
isoscalar interaction (xpp = Xpn) one has N n /N p = 1, and for a 
pure isovector interaction (xpp = ~ Xpn) one has N„/N p = — 1. 



III. WARD-TAKAHASHI IDENTITIES 

We first note the following commutation rela- 
tion [2(| between the angular momentum operators 
L ±x = 4= (L x ± iL v ) and a tensor operator T? k) of rank 
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k with spherical components q 
1 



.k 6 : 



r±l 

^ ' J (fe) 



-^v^fcTTF^iTyrg 1 . 



(mi) 



The commutator of Hq t (Eq. pi.4[ )) with the angular mo- 
mentum operators then becomes 



(III.2) 



and if we consider matrix elements of this identity be- 
tween non-interacting particle-hole states, we obtain the 
useful relation 



(a\L 



ill 



HQ? 1 !*) 



(III.3) 



which will be used in later Sections. (Here (at) € r.) 

Let us now discuss the Ward-Takahashi identity which 
follows from angular momentum conservation. We con- 
sider the time derivative of the 2-point function with ex- 
ternal Heisenberg operators Q\{t') and L x (t) = Lp{t) + 
L\(t). Using the Heisenberg equation of motion 



dL 



^=i[H,L l ]=iV3'£ Er Q 1 T , 



(III.4) 



and the equal time commutator [L 1 (t) , (t) ,] from 
we obtain the Ward-Takahashi identity 

j t {Q\T{Q?(t')L\t))\0) 

= iV3 J2(0\T {QV{t') Q{(t)) |0) e x - V3S(t - t')(Q° T ). 

(III.5) 

Let us define here the exact 2-point functions with 
one arbitrary operator (K — 1 component A 1 ) and the 
quadrupole operator (Q 1 ^): 

(0\T(Q?(t')Al(t))\0) = -iIl^(t'-t) 

= -i j due-^'^Ii^iuj). (III.6) 

Then the Fourier transform of the Ward-Takahashi iden- 
tity (|III.5|) can be expressed as 

^n? A L H = ]Tn? A %) £A -V3(Q° T ). (in.7) 



6 The definition used here for the spherical components of the 
angular momentum (£ = (L' x ± iL v )/^/2) differs in sign for 
the +1 component of any other vector (a 1 = —(a x + ia v )/\/2, 
a" 1 = (a x -iaV)/V2). Therefore, for the case T (1) = L in (JTTT7TJ, 
we have to use Tq} = T-^ 1 to get the correct commutation re- 
lation [L 1 ,^- 1 ] = L° = L z . 

7 We use the symbol H\ T for the exact correlators and the correla- 
tors in the chain (RPA) approximation, and ir\ T = S\ t tt t for the 
non-interacting ones. Note that ir^Q (u>) = 7t(oj) is the bubble 
graph of the previous Section. 



This is the basic identity which will be used in this pa- 
per. In the chain (RPA) approximation, the correlators 
(IIII.6|) can be expressed in terms of the reduced particle- 
hole t-matrix of Eq. ([II.13[) as follows (see Fig. 6): 

II^H = 5XT^ A (UJ) - JT r (u)Ux{u)^ A {u). (HI.8) 



For the case A = Q one can use Ea. (|TLT3)) to simplify 
this expression to 



■7r(w), 



1 + n(uj)x 

where we used the matrix notation of Eq. pi.14l) . 



(III.9) 



A. Identities for the Goldstone modes (uj — > first) 



In the limit uj -> (but finite e), Eq. (1111.71) leads to 
the following low energy theorem: 



(111.10) 



In the RPA, the correlator U qq (lu) is given by (|III.9[) . 
Inserting this form into pll.lOj) and multiplying from left 
by the matrix (1 + tt(0)x) we obtain 

m (P + x(Q°)) = (i + k(0)x) (Q°), 



where we used the self consistency relation pi.6[) to elim- 
inate e. (In this notation, j3 and (Q°) are considered as 
vectors in charge space.) We then obtain the identity 



tt t (0)/3 t = (Q° r ), 



(III.ll) 



which can also be shown directly by using the explicit 
form (|II.15|) of the bubble graph, see Appendix B. 

Using (jlll.llj) . the self consistency relation (|II.6|) can 

be rewritten as 



I3 T =e T -J2(XTXKx(0))f3 x 



(111.12) 



In the limit of exact rotational symmetry (e T = 0), this 
relation becomes 



Xpp7Tp(0) Xp„7T„(0) \(fip 

(0) J \ Pn 



(111.13) 



This equation leads to the condition Det (1 + x Tr (0)) — 
for a nontrivial solution. Comparing this with the pole 
equation (|II.16[) . we see that in the limit of exact rota- 
tional symmetry the self consistency relation guarantees 
the existence of a Goldstone pole (cjq — 0) in the K = 1 
channel. 

Let us now determine the vertex function for the Gold- 
stone modes. From pi.25[) and (|III.13[) we obtain the 
relation 



N P (0) _ I3 P 
N n (0) /3 n 



(III. 14) 
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FIG. 6: Graphical representation of the 2-point function Eq. ()III.8[l . 



Therefore the vertex function (|II.19[) for the K = 1 Gold- 
stone mode (n = 0) can be expressed as 

r™=°(x) = NQl(x)P T , (111.15) 

where (see (|III.14j) and (|II.21jl ) 

N T (0) = iV/3 T , ^ = X)/9?<(0). (HI.16) 

r 

The derivatives of the bubble graphs at ui = are ob- 
tained from pi.23|) and (|III.3|) as follows: 

a\LlH)\ 2 



(ai)Gr °" T (m)6T 



(111.17) 



Comparing this with the Inglis formula for the proton 
and neutron moments of inertia 12 711 



I T = 2 £ 



(111.18) 



we obtain the following important relations: 



#<(o) = v : 



AT 



3 



(111.19) 



where I = I p + I n is the total moment of inertia. The 
Goldstone vertex function (|III.15I) then takes the form 



Tr°(x) = JjQ 1 T {x)l3 r 



The vertex function for the K = — 1 Goldstone mode 
is obtained from (1III.20[) by replacing Q 1 — > Q^ 1 . It is 
easy to confirm that the sum or difference of the K = ±1 
vertex functions represents the change of the deformed 
mean field (U T (x) = Uq (x) — f3 T Q®(x) in the Hamiltonian 
(|II.4[1 ) under an infinitesimal rotation around the x or y 
axes |4j. 



B. Identities for the Q — L correlator (e — > first) 

For exact symmetry [e T = 0), the Ward-Takahashi 
identity T[) becomes 

c^n^H = -V3(Q° T ). (IH.21) 



In the limit uj 2 — > u> 2 , where u) n ^ is one of the 
nonzero solutions of the eigenvalue equation 16|) . the 
identity (|III.21I) gives 

lim {lo 2 - ujI) J2 n? A L M = o (w n + 0) . 

A 

(111.22) 



Inserting here the RPA form (|III.8[) and using the pole 
behavior of the reduced T-matrix pi,17[) . we obtain 

^7V A ( Wn )7T^ L K) = K ^ 0). (111.23) 



It is straight forward to check the validity of this relation 
by using the explicit form of the bubble graph ir® L (see 
Appendix B). 

Next, let us consider the limit to — > of (|III.21|) . For 
this purpose, we have to isolate the Goldstone pole on 
the l.h.s. Inserting the RPA form (IIII.8I) of the correlator 
UP L , and isolating the Goldstone pole by using pi. 171) . 
i.e., 



t\'\ 



N y (0)N x (0) 



+ (terms regular for uj — ¥ 0) , 

(111.24) 



the identity (|III.21|) in the limit uj — > becomes 



Iim^ ^7r T (0)iV T (0)iV A (0) 



(111.20) Using here pil.lip and pil.161) . we obtain 



lim, 



E 



fix = 



V3' 



V3{Q° T ). 
(11X25) 

(111.26) 



Again, it is easy to check this relation by using the ex- 
plicit form of the bubble graph n® L and the Inglis for- 
mula, see Appendix B. Actually, because both sides of 
pil.261) are one-loop quantities which consist of proton 
and neutron pieces, the identity (|III.26|) holds for sepa- 
rately for protons and neutrons: 



lim. 



(11X27) 



The identities (IIII.23[) and pil.271) will be useful for the 
discussion of transition matrix elements in the following 
Section. 
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\/2C0~n 




N T (a) n )Q^ 
->- 



CO, 



n 



For the case A = L = L p + L n , the identity (|III.23|) 
confirms that the total angular momentum operator can- 
not excite a state with finite excitation energy: 



{w n ,K = i\L 1 \o) = o K^o). 



(IV.6) 



FIG. 7: Graphical representation of the transition matrix el- 
ement, Eq. (|IV.l|) . 



IV. TRANSITION MATRIX ELEMENTS 

We first show that in the BS (RPA) framework the 
transition matrix element of the K = 1 component of any 
tensor operator A\ from the ground state to an excited 
state (excitation energy uj n and K = 1) is given by (see 
Fig. 7) 

(u> n ,K= 1|A^|0) = —^=W® A (uj n ) N T (w n ) . (IV.l) 

This formula reduces the determination of the summed 
transition strength 

B (A T ; K = 0^\K\ = 1) = 2 £ \(u n , K = 1\A\ |0)| 2 

n 

(IV.2) 

to a straight forward calculation of the Feynman diagram 
of Fig. 7. (The factor 2 comes from the contribution of 
K = -l.) 

To show pV.ip . we use the spectral representation of 
the exact correlator (jIII.6|> : 



= -E 



(0\Qy\n){n\A\\Q) (n\Q^ \0) (0\A{\n) 



ui — fl n + iS 



uj + Q n — i5 



(IV.3) 



Here 0„ are the exact excitation energies of the eigen- 
states \n) = \Cl n ; K = 1) of the Hamiltonian H. We then 
obtain for fl n > 

Jim^ (c 2 - fi 2 ) IP^V) - -2Q n <0|Qit|n)<n|4jO>. 

(IV.4) 



On the other hand, in the RPA we have from (|III.8I) and 
the pole behavior of the t-matrix pi,17[) 

hm ( w '-a£)n?>) 

UJ— >LJ n 

= -Tr T (uj n )N T (uj n )N x (uj n )TT^ A {oj n ). (IV.5) 



By comparing the r.h.s. of PV.4[) in the RPA (J7„ = oj n ) 
with the r.h.s. of (jIV.5l) . and noting that both expressions 
hold also for A — Q (where tt? = ir\), we immediately 
arrive at (|IV.1[) . 



In fact, the state L 1 10) has zero excitation energy because 
of [if, f 1 ] = 0, and pV.6[) confirms that this state is 
orthogonal to all states with finite excitation energy. 

For the transition matrix element pV.ll) of the oper- 
ator L T to the Goldstone mode, we can use the form of 
N T (0) from piI.19|) and the identity pil.271) to obtain 



(wo,Jif = l|4|0) = J^I T 



(uj -> 0) . (IV.7) 



If we sum over protons and neutrons, we obtain an iden- 
tity which follows also directly from angular momentum 
conservation, Eq. pll.26j) : 



(w 0> Jf = l|L 1 |0) 



up I 
2 



(wo-)-0). (IV.8) 



This relation shows that the state f ^O) is orthogonal to 
all RPA states, including the Goldstone mode. Actually, 
we will see in Sect. 6. A that Eq. pV.8l) is nothing but the 
normalization of the Goldstone state vectors. 



V. Ml SUM RULES 

As an application of the above formalism, we consider 
the inverse energy weighted (IEW) and energy weighted 
(EW) sum rules for the K = 1 component of the orbital 
magnetic moment (Ml) operator 8 



M 1 = 9lp + 9in Li, 



(v.i) 



where gi T are the orbital g-factors for r = p,n. (The free 
nucleon values are gf^ = 1, gf° Q — 0.) 

If we define the exact 2-point function with external 
Ml operators by 



n MM (w) = i / dre l " T (0|T(M lt (i')M 1 (i)) |0), 



= -2'£\(n\M 1 \ 



(V.2) 



where r = t' — t and we use the notations of Eq. pV.3p 
for the state vectors and energies, the IEW and EW sum 



Generally [3, a factor 1/3/(4-77) is included in the definition of 
the Ml operator. This factor is not included in our definition. 
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rules can be expressed as follows pjj 9 : 

^ s2E JWP = Bm n«w (v.3) 

n„>o 

S EW = 2 VKnlM^O)] 2 ^ = - lim w 2 n MM (w). 

' ^ UJ — hoc 

(V.4) 

To evaluate these sum rules in our BS (RPA) formal- 
ism, we introduce the correlator H^(uj), which has ex- 
ternal operators L^J and L\. The form of IT^(aj) in the 
RPA is (cf. Eq.lpLi])) 

U^(uj) = S tX ^ l (cj) - n^{Lo)t TX {uj)ir2 L {uj). (V.5) 

Here the first term on the r.h.s. is the non-interacting 
bubble graph. (For the explicit form, see Eq. (|A.6|) .) 

Concentrating first on the IEW sum rule, we note from 
the Inglis formula (|III.18[) that 7r^ L (0) is identical to the 
moment of inertia: 

*^(a, = 0)=2 V ™=/ T . (V.6) 

\oa) £r 

To evaluate the second term in (|V.5|) in the limit w — > 
0, we use the relation (|III.27|) . which shows that only 
the first (singular) term in the reduced T-matrix of 
Eq. (IIII.24[) contributes. Using the form of the normal- 
ization factors given in (|III.19[) , we obtain 

lim vr^H^HTrf'H = (V.7) 
The LL-correlator (IV.5I) for uj — > is then obtained as 

lim = 5 TX I T - (V.8) 

This result shows that the sum rule vanishes for the case 
where one of the external operators is the total angular 
momentum L, i.e., 

lmVlI#(w)=0. (V.9) 

UJ — 5-0 A ' 

A 

This is one of the many cases where the RPA-type corre- 
lations completely cancel the non-interacting (mean field) 
contribution, so as to satisfy the conservation laws (an- 
gular momentum conservation in the present case). 



9 The factor 2 in these expressions counts for the contribution from 
the K = — 1 component M~ 1 . (Equivalently, one can express the 
sum rules by M x or M y .) We also note that, at least in the RPA 
(see Eq. IIIV.7l l). the Goldstone term does not contribute to the 
spectral sum in I)V.2J| for finite uj, and therefore also not in the 
w -> limit of Eq. lfOt . 



Using (TVT8l) . we obtain for the IEW sum rule ((VT3]) 
Siew = lim V 5 £ T n^([j) 5 £A = —7-^ (gi.ivf , 

uj— >0 ' — » 7 

tA 

(V.10) 

where we define the isovector orbital g-factor by 

gtjv = 2 (#^> ~ Sfe) ■ ( v - n ) 

Before we continue to discuss the EW sum rule, we 
note the following two points: First, one can separate 
a term proportional to L 1 = + L n in the magnetic 
moment operator (|V.1[) according to 

M 1 = a L 1 + (gz >p - a) h\ + {g t>n - a) L\, (V.12) 

where a is any number. (The choice a = \ leads to 
the conventional separation into isoscalar and isovector 
parts.) From the result (|V.9|) it is clear that the first 
term in (|V.12j) does not contribute to the sum rule, and 
the contribution of the second term is independent of a, 
and given by (|V.10p . It is possible to choose a so that 
the RPA-type contributions vanish and the total result 
is given by the non-interacting correlator. This choice is 

a= I -fg tp + I -j gin = ge,is, (V.13) 

which leads to the following separation of the magnetic 
moment operator into "isoscalar" and "isovector" pieces: 

M 1 = g us i 1 + guv - ^) ■ (V.14) 

Both the IEW and EW sum rules discussed in this Sec- 
tion emerge exclusively from the second (isovector) part 
of Eq. (|V.14l) . We will see in Sect. 6 that this way to 
split the Ml operator follows naturally if one performs a 
minimal substitution in the effective Hamiltonians for the 
isoscalar (Goldstone) and isovector (scissors) rotational 
modes separately. 

Second, the IEW sum rule for the operator (2I n /I)Lp— 
(2ip//)Z/^ is often interpreted as the collective mass pa- 
rameter of the isovector rotatkm|12]. The result (|V. 10|) 
then confirms that this mass parameter is given by the 
"isovector moment of inertia" , which is defined as 

/iv = (V.15) 

Turning now to the EW sum rule (|V.4[) . we note that 
for the first term in the correlator (|V.5[) we have 

- lim uj 2 ^ l {uj)=2 V u ai \{a\Ll\i)\ 2 , 

= 3#7r r (0) = 3/MQ°), (V.16) 

where in the second equality we used the relation pil.3[) 
to express the result in terms of the bubble graph 7r T (0), 



and in the last equality we used the low energy theorem 
(IIII.11I) . To evaluate the second term in the RPA corre- 
lator (|V.5[) in the limit uj — > oo, we note that in this limit 
the t-matrix (III. 13|) becomes simply the 4-Fermi interac- 
tion constant X- Using also the form of the mixed bubble 
graph tt^(uj) given in (|A.5|) and the identity (|III.3|) . we 
obtain 

lim ^tt^H^aH^H 

= 3 (firltrW) XrX (&7r A (0)) = 3(Q°) X r A (Qj) ■ (V.17) 

Adding the pieces ()V.16[) and ()V.17[) we obtain 
- lim uj 2 U^(lo) 

= 3 ((Q° T )f3 T SrX + (Q° t )Xtx(QI)) ■ (V.18) 

Because of the self consistency relation (|II.6[) for exact 
rotational symmetry (e T — 0), we confirm that the ex- 
pression (|V.18j) vanishes if we sum over r or A, which is 
again a consequence of angular momentum conservation. 
In order to get the EW sum rule for the Ml operator, we 
can therefore discard the first term in (|V.14j) . By using 
again the relations (|II.6j) we finally obtain the following 
result: 

Sew = - lim u? V" g e T U^(ui)ge ,a, 

CU— ¥00 ' * 

T.X 

= -l2(g U v) 2 (Q° p )x P n(Ql). (V.19) 

This is essentially the result which has been obtained in 
Ref.[l4| by an explicit calculation of the corresponding 
double commutator. We also wish to mention that the 
EW sum rule for the operator {2I n /I)L x p ~ (2I p /I)Ll 
is often interpreted physically in terms of the restoring 
potential energy of the isovector rotation [12]. We will 
see in Sect. 6, however, that such an intuitive result for 
the restoring potential energy does not seem to emerge 
in our present BS (RPA) framework. 

Before we discuss the connection of these sum rules 
to observations, we make a comment on the orbital g- 
factors: In a more general approach, for example the 
Landau-Migdal theory [28j, the correlator (|V.2[) is rep- 
resented by the Feynman diagram in the first line of 
Fig. 8. The RPA-type correlations in the Landau-Migdal 
approach are included via the integral equation for the to- 
tal vertex, shown in the second line of Fig. 8. The driving 
term involves not the free but an effective magnetic mo- 
ment operator, which includes all processes (tensor cor- 
relations, meson exchange currents, etc) which are not 
taken into account by the RPA-type correlations. To 
a good approximation, this effective operator can again 
be represented in the state- independent form (IV.ll) . but 
now with effective orbital g-factors, which are different 
from the free ones. Therefore, in an RPA approach like 
our present one, the isovector g-factor in the sum rules 
(IV. 101) and (|V.19I) should be considered as an effective 
quantity, which is determined by data for magnetic mo- 
ments of neighboring odd-A nuclei. It is well known [29| 
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FIG. 8: Upper diagram: General representation of the corre- 
lator Eq. (|V.2[) . The black circle represents the effective irre- 
ducible particle-hole vertex, and the shaded square represents 
the full particle-hole vertex including the RPA-type correla- 
tions. Lower diagrams: Graphical representation of the BS 
equation for the full particle-hole vertex. The open square 
represents the irreducible particle-hole interaction. 



that this effective isovector g-factor gejv is larger than 
the free nucleon value (which is 0.5) by approximately 
10%. 

Let us now discuss the connection to observations: 
First, the 2 + state of the ground state rotational band is 
a classical example of an isoscalar type rotation. (Actu- 
ally, in Sect. 6 we will see how it emerges from the effec- 
tive Hamiltonian for the isoscalar rotational mode.) Its 
magnetic moment is determined only by the first term in 
(jVT4|) [30j . i.e., ju(2+) = gij S L with L — 2. Its g-factor 
is therefore given by 

g(2 + )^g LIS ^ I -fg p + I -fg n . (V.20) 

As a first estimate one can assume that Ip/I ~ Z/A and 
/„// ~ N/A, which gives 

Z N Z 

ge.is - -j gt, P + -rge,n = 7, (V.21) 

where the second equality has been derived rigorously 
from gauge invariance in a nuclear matter picture in 
Ref . [31j . We therefore obtain the familiar rotor value 

g{2 + ) = — . For example, the nuclei considered in the 

analysis of Ref. 1"> have Z/A ~ 0.4, and the measured 
g-factors are g(2 + ) ~ 0.3. To account for the difference, 
one has to take into account the effects of pairing, which 
enhance the neutron moment of inertia relative to the 
proton one [32], but qualitatively the rotor value is cor- 
rect. 

Second, it has been shown in Ref.[]j| that the IEW 
sum rule value (IV.10[) with I p /I ~ Z/A and I n /I ~ 
N/A, agrees well with the experimentally determined ra- 
tio B(Ml)/ui(Ml), where the transition matrix element 
B(M1) and the excitation energy oj(M1) refer to the ob- 
served low energy scissors mode. However, this agree- 
ment is only obtained if the isovector orbital g-factor 
gejv is assumed to have the same value as the isoscalar 
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one (ge.is)- As we discussed above, however, the isovec- 
tor orbital g-factor should be larger than the free nucleon 
value, while the isoscalar one is smaller. This seems to 
indicate an interesting problem which deserves further 
study. In particular, the effects of pairing 19] should be 
included in the present framework. Also, for a quanti- 
tative analysis one should extend the present model to 
include the effects of the nucleon spin, and investigate 
its role for the IEW sum rule, following for example the 
analysis of Ref.[14j for the EW sum rule. The effects 
of pairing and spin should be closely investigated before 
drawing conclusions on the problem of the IEW sum rule. 



VI. ROTATIONAL BANDS AND SCISSORS 
MODE 

In this Section we wish to derive the effective Hamilto- 
nians and the state vectors for the isoscalar and isovector 
rotational modes. Concerning the isoscalar mode, we can 
apply the formalism of Sects. 2 and 3 without further ap- 
proximations. For the isovector mode, however, it seems 
necessary to refer to the harmonic oscillator model in 
order to obtain analytic results. 

Let us first establish the connection between our for- 
malism, which is based on the BS equation, to the usual 
RPA formulation in terms of "forward" and "backward" 
amplitudes. This connection can easily be established if 
we return to Eq. pV.ll) for the transition matrix element, 
and insert the spectral form of the non-interacting bub- 
ble graph tt^ a , which is given by Eq. (jA.l[) for B = Q. 
In this way we obtain 10 



where O = Y. T °l with 

0\{u) n ,K = X) 

= [( Y r) m (ui n ) 4 a>i - (Z r ) m K) a] a a . 

(afi)GT 

(VI.5) 

Here and a are the single particle creation and anni- 
hilation operators. Second, one can derive an effective 
Hamiltonian for the RPA operators and O by follow- 
ing bonsonization methods based on the path integral [3| 
11 , i.e., introduce 0\ O as auxiliary quantities into the 
Hamiltonian and integrate over the Fermion Grassmann 
variables. The resulting Fermionic determinant can then 
be expanded in powers of and O. Since here we are in- 
terested only in the kinetic part of the effective rotational 
Hamiltonian (K — ±1 modes), we simply quote the re- 
sult, which is well known 0] 12 and can be motivated also 
by more intuitive arguments: 

ffeff = £0 + ^#eff(O, (VI.6) 

n 

where 

H eS (uj n ) = 

1 K=±l 

(VI.7) 



{w n ,K = l\A$.\Q) 



(az)£7 



(VI.1) 



Here we use the notation (O)- = (a\0\i) for the particle- 
hole matrix elements of an operator O, and the K = 1 
components of the RPA amplitudes are defined by 



(Y T ) m K) = 

( Z r) m (^n) = 



-1 N r (u n ) (Qr) m 
^UJn Ul n - UJ m + iS ' 
1 NrM (Q\) ia 

/2w„ u> n + U) m - i5 



(VI.2) 
(VI.3) 



We recall two results of the RPA: First, the states 
\ui n ,K = 1) can be expressed by 



\w n ,K = l) = Ol(u n ,K=l)\0), 



(VI.4) 



10 We follow the notations of Ref.p]], see in particular Eqs. (14.17), 
(14.26) and (14.32) of Ref.Q. 



A. Isoscalar rotational state 

Here we derive the form of the RPA amplitudes for the 
Goldstone modes (ujq — and K — ±1), and determine 
their contribution to the effective Hamiltonian (|VI.7[) . By 
using the form of N T (0) from piI.19[> and the relation 
pIL3]) , we find for the amplitudes Y and Z of (fVT2|) and 
(IVL3I : 



(Yr) m M 



(a\Ll\i) 



(z T ) m M = - 



i\Ll\a) 



(VI.8) 
(VI.9) 



11 Bosonization methods in the path integral formalism have been 
used in relativistic field theories to show the equivalence of 4- 
Fermi type interactions to Yukawa type interactions \3i , and also 
in nuclear structure physics to motivate the Interacting Boson 
Model [H. 

12 See Eq.(8.92) or (8.97) of Ref.Q. The constant Eq in iTVLo) 
corresponds to E^p — ^TrA in the notation of Ref. Q. 



11 



Inserting these forms into (|VI.5I) , we obtain the represen- 
tation of the Goldstonc state vector as follows: 



\lj ,K = 1) = O + (cj , K = 1)|0), where 



(VI.10) 



The creation operator for the second Goldstone mode 
(K = —1) is obtained by the replacement L 1 — > L^ 1 . 

We see that the creation operator for the Goldstone 
mode diverges as in the limit ujq —¥ 0. The nor- 

malization of the Goldstone state vectors, however, is 
finite, and obtained from (IIV.8[) and (|VI.10[) as 



(uo,K\uq,K) = 1. 

K=±l 



(VI.11) 



Inserting (|VI.10I) into the effective Hamiltonian (IVI.7|) , 
we obtain for the contribution of the Goldstone modes to 
the effective Hamiltonian: 



.ffeff(wo) = #rot(w ) 



{L x y + (L y ) 
21 



L' 2 



L 3 



21 



(VI. 12) 



Because this result was derived from the intrinsic Gold- 
stone modes, which are degenerate with the ground state 
of the spontaneously broken rotational symmetry, it cor- 
responds to the ground state rotational band. The point 
to note is that, while the Goldstone modes have zero in- 
trinsic excitation energy and correspond to the solution 
luq = of the BS equation, the operators 0^(ujq) and 
0(ujq) diverge as 1/y/uJo, and as a result the contribu- 
tion of the Goldstone modes to the effective Hamiltonian 
(|VI.7[) is finite and given by the collective rotational en- 
ergy (|VI.12[) . In a microscopic quantum theory of finite 
systems, the Goldstone modes are therefore by no means 
"spurious", but describe the rotation of the whole sys- 
tem around an axis perpendicular to the symmetry axis, 
with finite rotational energy. Only if one assumes the 
rotational part (|VI. 12[) from the outset, the Goldstone 
modes should be considered as "spurious" . 



B. Isovector rotational state (Scissors mode) 

Contrary to the isoscalar rotational modes discussed in 
the previous Subsection, it seems necessary to make more 
specific model assumptions in order to derive the corre- 
sponding expressions for the isovector rotational modes. 
Here we will assume the harmonic oscillator (h.o.) form 
Uo(x) — (Mcj 2 /2) r 2 in the mean field Hamiltonian 
(|II.4p . The sum Uo(x) — f3 T Q°(x) is then equivalent to 
a deformed harmonic oscillator potential. Although this 
model has already been used in similar contexts by many 
authors [6l. ll4ll36l . l37l ]. we discuss it in this Subsection and 
in Appendix C in some detail, because to our opinion it 
is highly interesting to see the isovector counterparts of 



the relations given in the previous Sections, even if those 
are more model dependent. 

In the h.o. model, the bubble graph of Eq. (|II.15p as- 
sumes a 2-pole form (see Eq. (jC.5[) ). which makes ana- 
lytic calculations possible. Here, in order to keep the 
equations as schematic as possible, we restrict ourselves 
to one h.o. shell (AJV = space), which corresponds 
to the first term in Eq. (|C.5[) . (The case of the full h.o. 
space can be found in Appendix C, and the main points 
will be summarized in the next Subsection.) In this ap- 
proximation, where all particle-hole states have the same 
excitation energy (eo T ), the bubble graph is given by the 
following one-pole form: 



7T r (w) = — 2- 



' uj 2 - e\ T + iS 
where e UT and Sq t are given by 



Sqt, 



eor = KPt, 

S r = J2\( a \QrW 



K 



-Or 



45 1 
16tT JdOi 

K 3 



(VI. 13) 



(VI. 14) 



G 



p T I T . (VI.15) 



The symbol in the sum (|VI.15j) indicates that only 
the AN = particle-hole states are included, and we 
used the identity (|III.3|) and the Inglis formula (|III.18|> 
to derive the second equality in (|VI.15|) . The low energy 
theorem (jlll.llj) relates Sq t to the quadrupole moment 
according to 13 



S0r = f (Q° T ) 



(VI.16) 



Using (|VI.14|) and (|VT.16]> . the self consistency relations 
(111.61) for the case of exact symmetry (e T — 0) can be 
expressed as follows: 



eop 
eon 



(VI. 17) 



Inserting the pole form (|VI.13[) into the eigenvalue equa- 
tion (|II.16[) . we get two solutions: The first one is the 
Goldstone solution (ojq = 0), and the second one is given 
by 



w l — e 0p { e 0p + ^XppSop) + eo n (eQ n + 2XnnSon) 



— —^Xpn (eOnS'op + GOpSon) 



(VI. 18) 



13 As is clear from Eq JC.51 1. the AN = sum So T is actually 
only half of the sum in the full h.o. space. This artifact of the 
restriction to one h.o. shell is formally remedied by replacing the 
quadrupole moment (Q?) by the core contribution (Q®) c , which 
is half of the total quadrupole moment, in all preceding relations 
of this paper. See Refs,([6l,l3Cll, l33lO for discussions on this point. 
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Using this solution, it is easy to calculate the proton and 
neutron normalization factors from Eqs. (|II.21l) . (1IL22[) . 
The result for the ratio, which corresponds to Eq. (jIII.141) 
for the isoscalar Goldstone mode, is 



N p (uji) u)\ 



-On 



'Op 



Son 



^OpXpp 



&0nXnp 

(VI.19) 

The important point to note is that, whereas for the 
Goldstone mode the ratio N n /N p is positive and close 
to +1, which corresponds to an isoscalar motion, the ra- 
tio (|VI.19|) is negative and close to -1, which corresponds 
to an isovector motion. For the individual normaliza- 
tion factors we obtain (cf. Eq. pil.19p for the Goldstone 
mode) 




(VI.20) 



(VI.21) 



Using these normalization factors and the expression 
(|VI.18[) for the excitation energy, it is easy to calculate 
the RPA amplitudes from (TVPj) and (jVL3i To illustrate 
the method, which is extended to the full h.o. space in 
Appendix C, we note that the operator 0\ of (|VI.5|) can 
be expressed as follows: 

Oi(. n ,K = l) = -^d 



1 



id 



Qr(0) 



(VI.22) 

where n = 0, 1. Here we used the identity pil.3[) . and 
defined the "low energy part" of an operator A as follows: 



A(0) 



o 



a] a 



(VI.23) 



In the AN = space, which is considered in this 
Subsection, we can identify these low energy operators 
with the full operators (A(0) = A). For the Goldstone 
mode (uiq — 0), Eq. (|VI.22p reproduces the general result 
(|VI. 10|1 . and for the oj-y mode we obtain 

|wi,iir = i) = ot( Wlj A- = i)|o) 

= (0{{u x ,K = l) + 0\ j {u l ,K=l))\Q), 

(VI.24) 

where the creation operator has been split into an 
angular momentum part ("L-part") and a quadrupole 
part ("Q-part") defined by 



0[(u u K = l) = 



1 



\/2uiIi V 



(VI.25) 




The isovector moment of inertia was defined in (|V.15I) . 
The L-part (|VI.25[) has the simple interpretation as the 
generator of an out of phase rotation of protons against 
neutrons, where the quantities 2I n /I and —2I p /I play 
the role of "effective charges" for protons and neutrons, 
which effectively remove the contribution of the overall 
in-phase rotation. The mode which is generated by this 
operator is therefore called properly the "scissors mode" . 
The Q-part (|VI.26|) . on the other hand, generates the 
quadrupole vibrations, and we will see below that it gives 
rise to the restoring force. (More general forms of those 
generators are given in Appendix C.) 

Including also the annihilation operators and the K = 
— 1 mode, we can summarize as follows: 



0(u u K=l) 
0{u!,K = -l) 



= 0[(lo 1 ,K = l) + Ol(u 1 ,K = 1), 

= O l (u 1 ,K = 1) + Oq(ui,K=1), 

= O l {u 1 ,K = 1)-Oq(ui,K=1), 

= Ol(u!,K = l)-Oi Q (ui,K = 1), 

(VI.27) 



where 0\(lj\,K — 1) and Oq(uji,K — 1) are given in 
dVL25l) and (IVL261) . and 



O l { Ui ,K=1) = 

Oq{ux,K=1) = 
1 



- ( 2I n _l 21p J- —1 

T L p T n i 



IV 



(VI.28) 



(VI.29) 



By adding the K = 1 and K = — 1 contributions to- 
gether, it is then easy to calculate the contribution of the 
uj\ mode to the effective Hamiltonian (|VI.7|) . We obtain 
the following result: 



(VI.30) 



where the rotational and quadrupole vibrational parts 
are given by 



(VI.26) 
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H tot (u i )=ui [<D[{u 1 ,K = 1)O l (ui,K = 1)+O l (uj 1 ,K = l)0\,(u u K = 1) 
Jfg(wi) = wi A" = 1)Oq(wi, A = 1) + Oq(uj u K = l)Oj,(wi, A = 1) 



(VI.31) 



(VI.32) 



By adding (|VI.12[I and (|VI.30|) . we obtain the total con- 
tribution of the isoscalar and isovector rotational states 
to the effective Hamiltonian: 

ffcff M + Heff(wi) = ffrot + HQ, (VI.33) 

where ffg = Hq{uj\) is given by (|VI32p . and # rot by the 
sum of (|VI.12[) and (1VL31I) : 

-ffrot = H To t(i0 ) + -ffrot(wi) 

_ (L'j? + {Llf | (L^ + (L^ (yl34) 

This is the kinetic part of the 2-rotor model Hamilto- 
nian [ll|. We therefore obtain the important result that 
the 2-rotor model is obtained from the RPA in a natural 
way by adding the effective Hamiltonians for the Gold- 
stone modes (isoscalar rotation) and the scissors modes 
(isovector rotation). 

We note that the two parts of the effective rota- 
tional Hamiltonian, given by the isoscalar Goldstone part 
(|VI.12I) and the isovector scissors part (|VI.31|1 . corre- 
spond exactly to the isoscalar and isovector parts of the 
Ml operator (|V.14[) . In the notation of first quanti- 
zation, this is seen most easily by making a minimal 
substitution, namely pi —¥ pi — gi p A{?i) for protons 
[i = 1....Z) and pj pj — gg n A(fj) for neutrons 
(j = 1,...N), in (|VI.12|) and (IVI.31|) separately Us- 
ing the form A(rk) — \ (b x ffcj, which corresponds 
to a constant external magnetic field, and expressing 
the magnetic interaction Hamiltonian in the form —B ■ 
M/(2M) = - {B X M X + Bv My) /(2M), gives the two 
terms in (|V.14p . Viewed in this way, the presence of 
the scissors part (|VI.31[) is necessary to give the correct 
coupling to an external magnetic field. 

The quadrupole vibrational part (|VI.32|) of the effec- 
tive Hamiltonian is positive definite (note that Xpn < 0), 
and represents the restoring force which acts against the 
proton-neutron oscillations in the present model. We 
note, however, that it cannot be reduced to a simple ge- 
ometric form, which is usually assumed in the 2-rotor 
model. 



C. Discussions 

In the previous Subsection, we have seen how the low- 
energy isovector scissors modes emerge in the simple ap- 
proximation of one major h.o. shell (A7V = 0). The case 
of the full h.o. space is discussed in detail in Appendix C, 
and we can summarize the results as follows: There are 
four solutions of the RPA equation, two at low energy cor- 
responding to the A7V = case discussed above, and two 
at high energy (A7V = 2). Each of these four modes can 
be represented similar to Eq. (|VI.24l) by an L-part and a 
Q-part, as shown in Eq. (|C.19p . However, the L-part now 
includes also the generator of quadrupole deformations, 
describing irrotational flow, in addition to the generator 
of rotations, as shown by Eqs. (|C.23|) and (IC.24I) . and the 
vibrational Q-part includes also the quadrupole operator 
in p-space in addition to the ordinary one in r-space, as 
shown by Eqs. (|C.25l) and (|C.26p . Concerning the L-part, 
for not too large deformations, the ordinary rotational 
term is dominant for the low-energy solutions, while the 
irrotational term is dominant for the high-energy solu- 
tions, as we explain in the last paragraph of Appendix 
C. 

Further information on the nature of these four collec- 
tive states can be obtained by considering their vertex 
functions and Ml and E2 transition matrix elements |14j. 
and one arrives at the following picture: One of the 
high-energy solutions (w2 in Appendix C) carries zero 
Ml strength, and is mainly of irrotational character. 
It can be identified as the K — 1 component of the 
isoscalar giant quadrupole resonance. The other high- 
energy solution (w 3 ) carries both Ml and E2 strength, 
and is also mainly of irrotational character. It is usually 
called the the high-energy scissors mode, or equivalently 
the K = 1 component of the isovector giant quadrupole 
resonance [l9j]. The nature of the low-energy solutions 
is essentially the same as we discussed in the previous 
Subsections, namely one is the isoscalar rotational mode 
(wo), and the other is the low-energy scissors mode (wi), 
which is mainly of rotational character. 

It is very interesting to note that this picture of ro- 
tational flow at low energy and irrotational flow at high 
energy is valid also for other many-body systems, like 
deformed metallic clusters [38| . deformed quantum dots 
[39| , crystals [4(J , and trapped Bose-Einstein condensates 
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[4lj . (The recent developments are reviewed in Ref.pjJ.) 
For example, in the case of deformed metallic clusters, 
the rotational flow corresponds to a rotation of electrons 
with respect to the jellium background, and the irrota- 
tional flow to a rotation of electrons within a rigid sur- 
face. In the case of trapped Bose-Einstein condensates, 
the collective modes are induced by an abrupt rotation of 
the deformed trap by a small angle, causing oscillations 
of the condensed atoms. This case is particularly in- 
teresting, because for normal (non-superfluid) gases one 
expects both the low-energy rotational and high-energy 
irrotational modes, while for the superfluid case one ex- 
pects only the high-energy irrotational mode because of 
the small moment of inertia. Experimental evidence that 
the low-energy rotational mode of trapped condensed 
gases indeed exists only above a critical temperature has 
been reported in Ref . [42j| . 



VII. SUMMARY AND OUTLOOK 

In this paper we used a simple field theory model based 
on a separable QQ interaction to gain analytic insights 
into the physics of rotational modes in deformed nuclei. 
Our essential tools were the Ward-Takahashi identities 
for angular momentum conservation, which we used to 
discuss the Goldstone modes associated with the spon- 
taneous breaking of rotational symmetry, in particular 
their vertex functions and decompositions into particle- 
hole components. In this simple model it was possible 
to derive analytically the ground state rotational band 
from the effective Hamiltonian for the isoscalar rotational 
modes. The isovector rotational (scissors) modes, on the 
other hand, correspond to a finite intrinsic excitation 
energy, and their properties depend on the mean field 
and the residual proton-neutron interaction. In order to 
obtain analytic results also for the isovector modes, we 
made use of the harmonic oscillator model for the spher- 
ical part of the mean field. It was then possible to derive 
the vertex functions, the decompositions into particle- 
hole components, and the effective Hamiltonian also for 
the scissors modes. By adding the effective Hamiltonians 
for the Goldstone and the scissors modes, we obtained the 
kinetic part of the 2-rotor model Hamiltonian, and also 
the potential energy (restoring force). 

An other important part of our analysis was the 
derivation of the inverse energy weighted and the en- 
ergy weighted Ml sum rules in the RPA, where we ob- 
tained analytic results without resorting to the harmonic 
oscillator approximation. We discussed those results in 
connection to recent experimental analysis on the scis- 
sors modes. We pointed out that other types of correla- 
tions, which are not included in the RPA (for example 
tensor correlations and meson exchange currents) should 
enhance the Ml sum rules, because those processes are 
known to enhance the isovector orbital g-factor. We 
pointed out an interesting problem in this connection: 
The experimental data for the inverse energy weighted 
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FIG. 9: Graphical representation of the bubble graph ir^ (w) 



sum rule seem to require that the isovector and isoscalar 
orbital g-factors are the same, while observations and 
theoretical calculations of magnetic moments clearly in- 
dicate that the isovector orbital g-factor should be larger 
than the isoscalar one. It is, however, necessary first to 
take into account the effects of pairing, and extend the 
model to include the spin of the nucleons, before one can 
arrive at firm conclusions. 

We finally mention that some parts of our analytic 
derivations should be possible for more general interac- 
tions, including spin-dependent and non-separable ones. 
The Ward-Takahashi identities should give an important 
guide for this purpose. 
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Appendix A: Form of bubble graphs 

In this Appendix we give the forms of various bubble 
graphs which appear in the main text. 

We represent a bubble graph 7nf A (w) with external op- 
erators B 1 ' and A 1 by Fig. 9. (We consider the case where 
these operators are the spherical K = 1 components of 
some tensor operators, since this is actually used in the 
main text.) 

Using the form pi. lit of the propagators and per- 
forming the integration over fco by residues, we get (see 
Eq. (III.15[) for the special case A = B = Q) 



= - E 



{a\A 1 \i){a\B 1 \i)* (i\A x \a){i\B x \a)* 



uj — uj a i + id 



uj + uj a i — iS 



(A.l) 



In order to combine these two terms, one changes the 
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sum over the single particle states (ai) in the second term 
to the time reversed states (ai), which have the opposite 
values of l z but the same energies for the axial symmetric 
case. Then one uses the property (see Eq.(A.22) of [l|) 



(i\0\a) = t (a\0\i) 



(A.2) 



where = +1 for a T-even operator (like the quadrupole 
operator in our case) , and to = — 1 for a T-odd operator 
(like the angular momentum operator in our case). As 
a result, the relative sign of the forward and backward 
terms is different for the cases where both operators have 
the same or the opposite T-symmetry, and one obtains 



7^(0,)= -2 ^ 



[a\A^)(a\B^\iY 



n 



- uSl; + iS ' 



(A.3) 



where f2 = u ai if tj± — is, and fl = uj if t& = —ts- 

For the cases needed in the main text, the bubble 
graphs are then obtained as follows: 

(A.4) 



1 



(ai) £t 



W 2 - W 2 ■ + 2(5 ' 



-2w ^ (ali^l^^alQ-^) 



1 



cj 2 - w 2 - + 



(A.5) 
(A.6) 



The two identical forms of the mixed bubble graph in 
(|A.5[) are obtained by expressing either the first or the 
second term in (|A.1|) by the time reversed states (ai). 
We also note that all bubble graphs which appear in this 
paper correspond to K = 1, although this is not indicated 
explicitly in our notations. 



Appendix B: Identities for bubble graphs 

In this Appendix we use the forms of the bubble graphs 
given in Appendix A to confirm various identities which 
are derived more generally in the main text. 



1. Eq. (|III.llj) 



We use Eq. (|A.4|) and the identity pil.3|) to write 
\{PH\Qi\0)\* 



7T T (0) = 2 £ 



(PH)Et 



UJph 



= 3gS £ »ph\(PH\I} 

' T (PH)er 



(B.l) 



where only in this and the following equation we denote 
the non-interacting (mean field) ground state simply by 
|0), the non-interacting particle- hole states by \PH), and 
their energies by ujph- Since (jB.ll) has the form of an 
energy weighted sum rule, it can be expressed as a double 
commutator in the standard way, and we obtain 



/WO) 



i 



[[H ,Ll],L- 1 ] |0), 
'<#,41|0) = <0|Q?|0), (B.2) 



where in the last two steps we used the commutation 
relations plTTTT) and (IITL21) . 



2. Eq. (|III.23|l 



If we insert the identity (|III.3[) into the first form (|A.5|) 
of tt® l and compare the result to the form (|A.4[) of n®® 
we obtain the following identity: 



Un? L (w) = V3/3 T (7T T (W) ~ 7T T (0)) . 



(B.3) 



Then, in order to show Eq. (|III.23p . we have to show that 
the following relation holds if ui = u> n is a solution of the 
eigenvalue equation (|II.16j) : 



(tt p (uj) - 7T„(0)) 



N n (u) p n 
N p (uj) f3 p 



(7r n (w)-7r n (0))=0. 



(B.4) 



Here the ratio of the normalization factors is given by 
(|II.22[) . and the ratio of deformation parameters by 
(MI. 141) . Using these relations, we can express (|B.4p 
solely by bubble graphs as follows: 



(lT p (uj) - 7T p (0)) 

(1 +Xpp 7r pM) 7r p(0) 



(7r„(w) - 7r„(0)) = 0. (B.5) 



(1 + Xnn^n(0)) Tt n (u)) 

More explicitly, the relation (|B.5[) can be written as 
Tr p (uj)ir n (Lu) - 7T p (0)7r n (0) 



Xpp 



7T p (a;)7r n (a;)7rp(0) + 7r p (w)7T n (u;)7r n (0) 
- 7Tp(cj)7rp(0)7r„(0) - 7r„(a;)7rp(0)7r„(0) 



= 0. 
(B.| 



In order to verify this relation, we use the eigenvalue 
equation (|II.16|) for u — ui n as well as for lu = 0. (Note 
that for cj = the RPA equation is equivalent to the self 
consistency relation, as we have shown in the main text.) 
This gives the following identity: 



(1 + XppTTpM) (1 + XnnK n (u)) _ TT p (u)TT n (Lo) 
(l+Xpp7Tp(0))(l + (0)) ^(0)tt„(0) ■ 



(B.7) 



It is readily seen that this relation is the same as (|B.6I) . 
This concludes the explicit verification of Eq. (|III.23I) . 
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3. Egs. pil.26|) and (|III.27|) 

Using Eq. (|A.5|) for the bubble graph tt^ l and the iden- 
tity (jTTTT3]) we obtain 



lim 



2 yv /; 1 



Using then the Inglis formula (|III.18I) we obtain pil.271) . 



This relation, which follows from the analytic forms of 
Sq t and S2 T given by Eq.(27) of Ref.[36], shows that the 
AiV = and AA = 2 excitations give the same contri- 
butions to 7r T (0). The low energy theorem then 
can be written in the form 



S r = f (Q° T ) 



(C.8) 



Using (|C3[) and (|C.8|) . the self consistency relations (|II.6[> 
take the form 



Appendix C: Deformed harmonic oscillator 

In this Appendix we review some formulas @, [ToL ITil 
HH, H3] for the RPA calculation with deformed h.o. mean 
fields. 

Using U (x) = Mw 2 r 2 /2 in the mean field Hamiltonian 
(|II.4p . the sum Uo(x) ~ (3 T Q°(x) becomes a deformed h.o. 
potential with frequencies 



^z,r = - -<5 T 



(C.l) 



e 0p — ^4x pp S'op — 4x P n Son, 
eon = — 4x„„ 5on — 4x„ p 5q p . 



(C.9) 
(CIO) 



Inserting the form (|C.5I) of the bubble graph into the 
eigenvalue equation (|II.16|) . we can calculate the col- 
lective excitation energies in this model. Besides the 
Goldstone solution (luq — 0), there are three solutions 
(wi, 012,^3) with positive energy, which are determined 
by the following cubic equation[6J in x 



,2. 



with the coefficients 



bx — c = 0, 



(C.ll) 



where the dimensionless deformation parameters 8 T are 
related to the (3 T of the main text by 



5 T = —Pt 

to 



(C.2) 



and K is defined in Eq. (|VI.14p . In this simple model, two 
types of particle-hole excitations contribute to the bubble 
graph of Eq. (|II.15|i . corresponding to excitations within 
one h.o. shell (AA~ = 0) and across two shells (A A" = 2). 
The corresponding excitation energies are given by 



G(] T — LJ XtT — ^z,T — &t ^ ■ 



62r = ^x,t + W Z , T — 2w ( 1 - -5 T 



The bubble graph 151) then takes the form 

7T T (w) = — 2 ( ° T „ S T H ? 2r 2 for ) ! 



(C.3) 
(C.4) 



■ : Or 



'2t 



2 0t , e 2r A ^Or 



9 2 9 9 

w z - eg T - e^ T / e T 



(C.5) 



Here the quantities 

s 0t = ^IHq^ 



s 2t = ]T|HQ^)| 2 , (c, 



denote the sums over the AA" = and AA" = 2 particle- 
hole states, and in the second equality of (IC.5|) we used 
the relational 



Sot _ S 2t 



(C.7) 



a=(e 2 0p + el p )W p (2) + (j)^n), 
b = e 2 0p e%W p (4) + (p -> n) 



(C.12) 



+ J (4 P + 4 P ) {4n + tin) (1 + W P (A) + W„(4)) , 

(C.13) 



4 
1 



c = 2 e o P e 2 P (eon + eln) Wp( 4 ) + (p -> »)• (C14) 



Here we defined (for fc = 2, 4) 

Wr(jfc) = l + fc^SW- 

eor 

For the case where the proton and neutron deformations 
can be assumed to be equal {8 p = 5 n , which implies eo p = 
6Q n = eo and e2 P = e2n = ^2), simple analytic solutions 
of (|C.11|) exist: One can be obtained by noting that for 



lo 2 



el 



(C.15) 



the bubble graph (|C.5[) has the same value as for ui = 0. 
Since wo = is a solution of the eigenvalue equation 
because of the self consistence relations, (|C.15|) is also 
a solution. To further understand the physical nature 
of this solution, we note that because of ir(uJ2) = ""(0), 
the identity (|B.3|) shows that TT QL {bj 2 ) = , i.e., the Ml 
transition matrix element (B(M1) of (I1V.1[) ') vanishes for 
this mode. On the other hand, B(E2) is non-zero .14]. 
and this mode can therefore be identified as the K = 1 
component of the isoscalar giant quadrupole resonance 
[14| . The remaining two solutions can then be found 
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by solving simple quadratic equations, 
lowest order in 5: 



wi = e \ 1 



2 + 6' 



w 3 = e 2 \/l + -, 



One obtains to 

(C.16) 
(C.17) 



where 



Xpp - Xpn _ X( T = jj 



(C.18) 



is the ratio of the isovector to the isoscalar interaction 
strength (b > 0). Analytic solutions can be worked out 
also for the case of different proton and neutron deforma- 
tion parameters, although the expressions become quite 
long. To summarize, there are four solutions of the RPA 
equation, where ujq, uj± are the "low-energy" solutions, 
and u>2, W3 the "high-energy" solutions. For each solu- 
tion, one can determine the normalization factors N p (u n ) 
and N n (uj n ) from pL2T1) and pL22]) . 

Let us outline here the calculation of the creation oper- 
ator (|VI.5[) for any of these modes, following the method 
explained in Refs.fiol. l36j 14 : If we add the contribution 
of the AN = 2 excitations to Eq. (|VI.22p of the main text, 
we obtain 



Q\w n ,K = l) 

E 



N{w n ) 



0\{u n ,K=l) + 0Uu n ,K = \). 



(C19) 

For simplicity, we write the following expressions in first quanti- 
zation and omit the distinction between protons and neutrons. 



Here ujq, ...o>3 denote the RPA eigenvalues, and in ad- 
dition to the "low energy part" of an operator A, which 
was defined in (IVI.23[) of the main text, we also define 
the "high energy part" as 



A(2)=J2[(A) ai 



(C.20) 



In the harmonic oscillator model, explicit forms of the 
operators L 1 (m) and Q (m), where m — 0,2, can be 
derived as follows: We have 



Q 1 = - 



[(x + iy)j p zj - (p x + ip y ) J 




(C.21) 



(C.22) 



where j labels the nucleons. If we express (|C21|) and 
(|C.22[) in terms of the standard creation and annihila- 
tion operators a\ and (k = x,y,z) for each parti- 
cle, we obtain two kinds of terms: The first kind in- 
volves products a\a,i and a\a z with i = x,y, and the 
second kind involved products a\a\ and a z a,. It is clear 
that the first kind of operators contributes exclusively 
to AN = excitations (operators A(Q)), and the sec- 
ond one exclusively to AN = 2 excitations (operators 
A(2)). Then, for each part A(m) separately, one can re- 
expresses the creation and annihilation operators by the 
original position and momentum operators. In this way 
one obtains the decompositions L 1 — i 1 (0) + L 1 (2) and 
Q 1 = Q\0) + Q 1 {2), where 



(d> x + Cj z ) J ! . UJ, ■ ^ 



UJ X + UI 



(u> x + uj z f 
4£> X L) Z 



1 LO x - 0J Z 

(1) U) x + ij z 



(1) (1) 
T j X Pj 



(2) 



(C.23) 



L 1 (2) 



(uj x - LJ Z ) 2 J j _ .U!, r - ^ 



ACj x lo z 



(u> x - co z y 



(1) U x - UJ Z 



(1) (1) 
r ) x Pj 



(2) 



(C.24) 
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and 



> 1 (0) = ^Q 1 - 



8tt 2M 2 



P*j (Px +iPy) i 



1 



(2) M 2 U) x &z 



1 



(2) M 2 G) X U) Z 



(1) (1) 
P • x V ) 



(1) (i) 
Pj x Pj 



(2) 



(C.25) 



(2) 



(C.26) 



Here la™ x 

^^(fe) denotes the tensor product, with 

rank k and spherical component q, of two vectors a and b, 
according to the definitions of Ref. (2(| . The forms given 
above exhaust all possible one-particle tensor operators 
which can be formed from the position and momentum 
operators. 

Inserting these forms into (IC.19I) . we arrive at the fi- 
nal expression for the creation operator of each mode. 
We also note that the relations (j VI. 2 71) of the main text 
are still valid with the above extended operators, and 
therefore also the expressions given in the first lines of 
Eq. (rVT31|) and (TVT32j) remain valid. 

We see that, in addition to the two opera- 
tors L 1 and Q 1 , the presence of the high energy 
modes leads to two more types of operators, namely 



[ r (i) x p W 
operator U 



1(2) 
.(1) 



and 

x p^ 1 



[pW x p 



(Dl 



1(2)- 



In particular, the 



] , 2 « is the generator of quadrupole 



deformations [1C|, and describes irrotational flow. This 
is easily seen by noting that, for example, the part 
in { . . . } in the first line of Eq. (|C.23j) , which corre- 
sponds to the motion around the x axis, is given by 

^ (VP* - ZPy) } + - {VP, + ZPy) j 



It generates the 



displacement oc e x x r H — Vj/z I of the volume element 



of the liquid, which consists of a rotational and irrota- 
tional part. We also note that, in contrast to the angular 
momentum, the generator of quadrupole deformations is 
not a symmetry transformation of the Hamiltonian. 



For not too large deformation, we see from (|C.23|) and 
(|C24|) that the rotational term is dominant in the low- 
energy part -^(O), and the irrotational term is dom- 
inant in the high-energy part i 1 (2). Going back to 
Eq. dC.19p . this implies that the low enery solutions cor- 
respond mainly to rotational flow, and the high energy 
solutions mainly to irrotational flow. Further informa- 
tion on the character of these modes is obtained by con- 
sidering their vertex functions and Ml and E2 transition 
matrix elements 14], which leads to the following picture: 
The wo mode is the isoscalar rotational (intrinsic Gold- 
stone) mode, and the 0J2 mode is the K = 1 component 
of the isoscalar giant quadrupole resonance, as noted al- 
ready above. The remaining modes u>± and W3 are of 
isovector type, and carry both Ml and E2 strength. The 
loi mode is mainly of rotational nature and is called the 
low-energy scissors mode, while the W3 mode is mainly 
of irrotational nature and is called the high-energy scis- 
sors mode or, equivalently, the K = 1 component of the 
isovector giant quadrupole resonance (l9j. 
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